Skyrmion Lattice in Two-Dimensional Chiral Magnet 
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We develop a theory of the magnetic field-induced formation of Skyrmion crystal state in chiral 
magnets in two spatial dimensions, motivated by the recent discovery of the Skyrmionic phase of 
magnetization in thin film of Feo.5Coo.5Si and in the A-phase of MnSi. Ginzburg-Landau functional 
of the chiral magnet re-written in the CP^ representation is shown to be a convenient framework 
for the analysis of the Skyrmion states. Phase diagram of the model at zero temperature gives a 
sequence of ground states, helical spin Skyrme crystal — > ferromagnet, as the external field B 
increases, in good accord with the thin-film experiment. In close analogy with Abrikosov's derivation 
of the vortex lattice solution in type-Il superconductor, the CP^ mean-field equation is solved and 
shown to reproduce the Skyrmion crystal state. 

PACS numbers: 



I. INTRODUCTION 

Skyrmions, originally proposed as a model for baryons 
in nuclear physics^, were first realized experimentally 
in the condensed-matter system of quantum Hall ferro- 
magnets near integer filling factor v « 1^'^. The elec- 
trons spontaneously form a fully polarized ferromagnet 
a,i V — 1 due to the exchange interaction, while slightly 
away from it the spins of extra (or lack thereof) electrons 
organize themselves into an intricate Skyrmionic struc- 
ture as a result of the competitive interplay between the 
Zeeman and the Coulomb interactions^. It was further 
suggested that these Skyrmions may condense into a crys- 
talline form**. Skyrmions as quasiparticle excitations in 
quantum Hall ferromagnets are by now well established 
experimentally'^. Recent experiments also find support 
for the crystallization of Skyrmions in the quantum Hall 
system^ . 

More recently, a strong case for the formation of 
Skyrme crystal (SkX) state was presented in the A- 
phase of a metallic ferromagnet MnSi^ and other com- 
pounds of B20 structure^'*. Here, the main evidence 
comes from small angle neutron diffraction data which 
exhibits a clear hexagonal pattern consistent with the 
triangular lattice arrangement of Skyrmions. Subsequent 
observation of the anomalous Hall effect in the A-phase 
gives additional support to the Skyrmion lattice picture^ . 
On the other hand, Monte Carlo simulation of the two- 
dimensional classical spin model with the Dzyaloshinskii- 
Moriya (DM) interaction and magnetic anisotropy pre- 
dicted the existence of various Skyrmion crystal phases 
as the ground state^*^. Just after this theoretical pro- 
posal, a real-space observation of the Skyrmions form- 
ing a crystalline phase has been provided by the Lorentz 
TEM imaging technique in the thin film of a non- 



centrosymmetric magnetic crystal Feo.5Coo.5Si under a 
perpendicular magnetic field^^. By now there is accu- 
mulating experimental evidence stating that Skyrmion 
lattice is a natural occurrence in magnetic metals lack- 
ing inversion symmetry. For such crystals the role of 
the DM interaction is to convert a k = ferromagnetic 
ground state to a helical phase at a finite ordering vector, 
k 7^ 0, with the magnitude |k| fixed by the ratio of the 
DM exchange over the Heisenberg exchange energies. In- 
fluence of thermal fluctuation and/or the magnetic held 
can favor the formation of multiple spiral phase over a 
single spiral (In this paper we use the words "spiral" and 
"helical" interchangeably). When imaged in real space, 
the multiple spiral phase is none other than the Skyrmion 
lattice^^i"'". 

We should emphasize the contrasting behavior of the 
chiral magnet realized in two (2D) and three (3D) spatial 
dimensions. For 3D systems such as examined in Refs. 
6-8, the Skyrmion crystal forms over a small window of 
flnite temperature just below the paramagnetic transition 
known as the A-phase. On the other hand, for quasi-2D 
chiral magnets recently synthesized with the thickness 
smaller than a single spiral period"'^^, the Skyrmion phase 
occurs over a much larger temperature range extending 
nearly up to the paramagnetic transition and down to 
the lowest temperature measured, T « 5K, strongly sug- 
gesting that the Skyrme crystal is the ground state. Free 
energy analysis of the 3D model indeed shows stability of 
the Skyrmion phase at a flnite temperature, but not at 
zero temperature^"^. As shown in this paper, the 2D case 
supports the Skyrmion phase even at zero temperature. 

Bogdanov and collaborators had earlier investigated 
the instability of the helical spin state to the spontaneous 
creation of Skyrmionic spin texture due to magnetic field 
in chiral magnets^^'^"^. Although not explicitly empha- 
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sized at that time, present experimental situation^^^'^^ 
shows that the dimensionahty plays a great role in sta- 
bilizing Skyrme crystal phase at low temperature. Our 
paper therefore re- visits the issue of the phase diagram of 
the chiral magnet under the influence of magnetic field 
first pioneered by Bogdanov et al, in light of these re- 
cent developments. We base our calculation explicitly 
on two dimensions, assuming spin models with zero spin 
anisotropy such as known to be the case for Fco.5Coo.5Si. 
As a theoretical advance, we show that the Ginzburg- 
Landau action for the chiral magnet can be cast in a 
simple form using the CP"'^ representation of the classi- 
cal spin. The advantage offered by the new representa- 
tion becomes evident when we consider the lattice case, 
where it becomes possible to derive the solution for the 
Skyrmion lattice state by exploiting the analogy to an- 
other celebrated example of the lattice of topological de- 
fects - Abrikosov vortex lattice. 

This paper is organized as follows. In Sec. II we in- 
troduce CP^ formulation of the Ginzburg-Landau (GL) 
energy functional for a chiral magnet. In Sec. Ill we de- 
rive the zero-temperature phase diagram of the 2D chiral 
GL model using two complementary approaches. Both 
lead to the ground state evolution with the external field 
in excellent accord with the recent experiment. Ana- 
lytical solution of the Skyrme lattice state is presented. 
Finally we conclude with a brief, qualitative discussion 
of thermal fluctuation effects for the Skyrmion crystal 
and summarize in Sec. IV. One can read this paper as 
a companion article to Ref. 11, where the experimental 
findings and Monte Carlo results were reported. 

II. CP^ FORMULATION OF CHIRAL 
GINZBURG-LANDAU THEORY 

We are concerned in this paper with magnetic sys- 
tems lacking an inversion symmetry such as MnSi and 
Fei_a;Coa;Si, where the magnetic interactions arc char- 
acterized by a ferromagnetic exchange J and a weaker, 
Dzyaloshinskii-Moriya coupling D. Hereafter, we take 
the length of a structural unit cell as unity. The competi- 
tion of the two interactions leads to a helical spin ground 
state with the pitch vector of length k = D/J. The effec- 
tive continuum theory for such chiral magnet is provided 
by the Ginzburg-Landau (GL) energy functional^^"-^^ 

J^[n] = ^ ^(5^n) • {d^n)+Dn ■ V x n-B • n. (1) 

n 

The non-linear cr-model given in the first part describes 
the excitation of a conventional ferromagnet while the 
second term gives the effect due to the DM exchange. As 
we arc interested in the effectively two-dimensional chiral 
magnet such as recently synthesized with Feo.sCoo.sSi^^, 
the spatial derivatives run over the two-dimensional 
plane, /i = x,y. The total energy is the spatial in- 
tegral of J^: F[n] = J (Pr J^[n]. It is also known 



that spin anisotropy does not play an important role in 
Feo.5Coo.5Si, which is therefore omitted from the theory. 

According to the Monte Carlo simulation^" and also 
the recent experiments^, the spiral spin states at B = 
is supplanted by the hexagonal packing of Skyrmions 
when B oriented perpendicular to the thin film plane 
exceeds a certain strength, Bd- At still higher field 
strength B > Bc2 the Skyrme crystal state gives way to 
fully polarized ferromagnet. Both transitions are first- 
order, as evidenced by the presence of the co-existence 
region around both B^ and Bc2^^- For the interme- 
diate phase Bel < B < Bc2 the triangular lattice of 
Skyrmions identified by the Lorentz TEM bears strik- 
ing resemblance to the Abrikosov vortex array in type- 
II superconductors^^'S^. In this paper we correctly re- 
produce the observed first-order transitions of magnetic 
phases under increasing magnetic field, by comparing the 
energies of the three candidate states: H (helical), SkX 
(Skyrme crystal), and FM (ferromagnetic). 

Previous theoretical approaches*'' ^^^^^ were based on 
the 0(3) representation of the spin vector, n = 
(sin^cos(?i, sin sine/), cos ^). In this work, we adopt in- 
stead a "complex" description in which n is replaced by 
a pair of complex fields 

^-[z2j~[ sin(e/2) j' 

obeying the unit-modulus constraint z^z = 1^^. Map- 
ping to the spin vector is through n = z'^crz. A well- 
known identity allows the re-writing of the nonlinear 
sigma model in the CP^ language (no sum on /u)^^ 

lid.n) ■ {d,n) ^ (d^z^) ■ (d,z) - Al, 

A, = -^[zHd,z)-id,z^)z]. (3) 

The Skyrmion density is related to the magnetic field 
associated with the vector potential by ^n- {d^n xdyii) = 

dxAy — dyAx- The chiral, DM term in the energy (1) 
under the CP^ mapping becomes 

n- (Vxn) = -2n- A-iz\a ■ V)z+i{Vz^) ■ crz. (4) 

When Eqs. (3) and (4) are combined, the energy density 
re-written in the CP^ representation takes on a succinct 
form 

7-[z] = 2J^(^^zy(D^z) -B-ztcrz, 

= di^-iA^-zKa^. (5) 

The derivative Z?^ is a 2x2 matrix due to the non-zero 
DM exchange, k = D/2J. This CP^ version of the energy 
density of a chiral magnet and the saddle-point equation 
which follows from it form the basis of the subsequent 
analysis^^. 
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To orient the readers, we first discuss how to write 
down a conventional magnetic state in the CP^ language. 
The ferromagnetic state for instance is written zpM = zq, 
an arbitrary constant spinor. The helical spin state - 
the ground state of Eq. (1) at zero magnetic field - is 
obtained from FM by a position(r)-dependent twist of 
the spin orientation, 



ZH 



Zq. 



(6) 



The Pauli matrix cr appears above. It can be chocked 
that the associated spin configuration nn = zJjCtzh is 
indeed orthogonal to the propagation vector direction 
k = k/|k| and rotates with the pitch k = 2k — D/J, 
provided the initial spin orientation no = zJctzq is per- 
pendicular to k. The vector potential A for the helical 
spin state is identically zero. 

An isolated Skyrmion in two dimensions with the spin 
pointing down at the origin r = and up far away, r — )• 
oo, is given the CP^ expression^^ 



state if wc consider the full free energy in Eq. (1). Al- 
though the exact spin configuration will be somewhat dif- 
ferent from that given in Eq. (8), the topological prop- 
erties remain the same. Furthermore, the radius R is 
uniquely determined by the ratio of DM interaction and 
the exchange energy as shown in the; next section. For 
Skyrmions in the quantum Hall system the radius was 
fixed by the relative strengths of Zeeman and Coulomb 
interactions^. 
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FIG. 1: (a) A typical (anti-) Skyrmion configuration given by 
Eq. (8) witli R = 2. (b) zi component of the CP^ Skyrmion 
zsk- When expressed as a planar spin (Re[2i], Im[«i]), it is an 
anti-vortex. 



zsk = e*(«W/2)(--fxi)^^^ 



(7) 



where zq = 



and f = (cos sin 0) is the ra- 



dial vector (Wc use (6 for the azimutlial angle of the 
spin vector n, and (p for the azimuthal angle of the co- 
ordinate r = {x,y): tarup = y/x). The polar angle 
9{r) is a function of the radial coordinate r such that 
0{O) = TT, 9{(X)) — 0, and smoothly varying in between. 
Vector potential for the single- Skyrmion configuration is 
Ask = i'P/''') sin^(6'/2). As a concrete example of a single 
(anti-) Skyrmion one can take the following case, 



r'^+R^' 



n +in 



V 



2iR 



x+iy 
r^+R'^' 



(8) 



and the associated CP^ expression zgk = 

^j^^i I _^ 1 • Note that the zi part of the 

CP^ Skyrmion solution behaves exactly as would a 
(anti-)vortex: depleted to zero at the core, winding 
about the origin as ~ e""'', and reaching a constant 
unit magnitude when r/R ^ 1 (Sec Fig. 1). At least 
pictorially it is clear that U(l) vortex configuration in 
the CP^ representation (in zi for instance) implies a 
Skyrmionic spin configuration in n through the mapping 
n = z^crz. In the next section, we will show that the 
mapping from a U(l) vortex to the real-spin Skyrmion 
can be extended to the lattice case as well. 

The spin configuration Eq. (8) is the saddle point solu- 
tion of the first, Hcisenbcrg-only term of the free energy 
given in Eq. (1). It is an excited state, with the energy 
E = An J independent of the radius of the Skyrmion R, 
and carrying a non-trivial topological number""^^. On the 
other hand, the Skyrmion state can become the ground 



III. PHASE DIAGRAM 

In this section we address the phase diagram of the 
model, Eq. (1) or Eq. (5), as the external field B = 
Bz, always assumed to be in the positive ^-direction 
(B > 0), increases from zero. Experimentally, the low- 
temperature phase evolves as helical spin, SkX, and FM 
with increasing field strength. We adopt two approaches 
for calculation of the energies of the respective phases, 
both of which successfully reproduce the observed phase 
evolution. In addition, we obtain analytic results for the 
two critical fields Bd and Bc2 {Bd < Bc2), each refer- 
ring to the first-order critical field separating H from SkX 
(Bel), and SkX from FM (Bc2). 



A. Vciriational analysis of a single Skyrmion 

In this subsection, we regard the Skyrmion lattice as 

the close-packing of individual Skyrmions of radius R 
forming a triangular lattice. The local spin orientation 
{9, (j)) of a single Skyrmion depends on the local coordi- 
nate (r, (f) as (f> = cp — tt/2 and 9 = 9{r). The total free 
energy of a single Skyrmion reads 



2J 



2Trrdr 
1 



ld9_ 
2dr 



■ sm tl cos ( 



-h ^sin^6'-/3(cos6i- 1) 
4r^ 



(9) 



where (3 = B/{2J), and the FM state is chosen to have 

the free energy zero. Ferromagnetic configuration is en- 
forced in the outermost region by the upward magnetic 
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FIG. 2: (color online) Dependence of the free energy density 
-Fskx on the radical cutoff R at magnetic field /? = I (red 
curve) and /3 = (black curve), respectively. We choose k — 
1. The angle measured from the «-axis as a fuction of the 
distance r from the center of the Skyrmion for two different 
R's are shown in the insets. Optimal _Ro for P — 1 leads to 
an almost linear dependence of 9{r) on r, while larger R = B 
leads to a long ferromagnetic tail outside the Skyrmion core 
region with radius = Rq. 



typical Skyrmion size as well as half of the optimal inter- 
Skyrmion distance, and the Skyrmions are close-packed. 
Applying the dimensional analysis to Eq. (9), one finds 
r has the same dimension as 1/k ~ J/D. Therefore it is 
expected that Rq would be proportional to which is 
confirmed numerically as well. 

When the cutoff R exceeds the optimal distance Rq, 
the free energy starts to increase from the negative min- 
imal value and approaches zero as i? — > oo. From the 
corresponding configuration shown in the right inset of 
Fig. 2 one can see that a.t R > Rq, the variation of 
9{r) is qualitatively different from that a,t R — Rq. A 
ferromagnetic tail with 9 = appears over the region 
between Rq and R. In usual conventions, only the core 
region where the variation in 9{r) is nonzero is referred 
as a Skyrmion. Therefore two neighboring Skyrmions are 
well-separated by the intervening FM phase in this case, 
and the Skyrmions are not close-packed. 

The other competing phase observed besides the 
Skyrmion lattice is the spiral configuration at small mag- 
netic field. This helical phase is the exact ground state 
of the Hamiltonian containing the DM interaction when 
the external field is absent. To see how the free energy of 
the helical state varies with field, consider the spiral spin 
propagating along y direction with (j) = and 9 = 9{y) is 
a function of y. The free energy for half period of spiral 



field, so that 9{po) = 0. On the contrary, 9{Q) = tt due 
to the geometric nature of Skyrmion. In numerical cal- 
culations, a hard radical cutoff R is introduced such that 
9{r) = for r > i?. In physical terms R is half the inter- 
Skyrmion distance in the Skymre crystal. The Skyrmion 
lattice observed experimentally can be constructed as the 
close-packing of non-overlapping, individual Skyrmions 
in a trianglular lattice. One can write down the total 
free energy of the Skymre crystal state as 



kX 



(10) 



L being the sample size. For each cutoff i?, we can ap- 
ply the numerical variation to minimize the free energy 
Fsk for a single Skyrmion. The equilibrium configura- 
tion of the whole lattice should minimize i^skXj hence 
the free energy density functional F's,\i[9{r)\/ R^ by opti- 
mizing 9{r). 

Figure 2 shows dependence of the free energy i^skx on 
the cutoff i?, where the sample size is normalized. It 
clearly shows that when R increases from zero, -Fskx de- 
creases dramatically, and reaches a minimal value at Rq. 
This point is exactly the ground state of the helical mag- 
net at zero temperature where the compromise between 
Zeeman, Heisenberg and DM interactions are reached to 
the maximum degree. The configuration of 9{r) is shown 
in the left inset of Fig. 2. One can find that 9(r) varies 
almost linearly in the whole region from to Rq, and 
reaches zero at Rq. In this case, Rq defines both the 



2JL 



I 



dy 



1 d9 



-^(cos6'-l) 



(11) 

with the boundary condition 0(0) = tt and 9{i) = 0, 
where 2(. is the period of a spiral. Similar to the Skyrmion 
lattice case, one can derive the optimal I by minimizing 
this free energy functional to get the configuration for a 
spiral lattice. 

The free energies of the three phases (helical, SkX, FM) 
for K = 1 obtained by variation are shown in Fig. 3. It is 
explicitly shown that the spiral state has the lowest en- 
ergy among the three configurations when the external 
magnetic field is small. In the large field limit, ferro- 
magnetism is energetically favored. In the intermediate 
region, one finds the Skyrmion phase as a compromise be- 
tween the DM and Zeeman energies. This is because the 
Skyrmion configuration keeps the spiral structure inside 
and ferromagnetic state outside, and hence can gain both 
the DM energy and magnetic field energy. This result is 
quite consistent with the Monte Carlo simulation^" and 
the experiments^. The lower (/3ci) and upper (/3c2) crit- 
ical field strengths separating the Skyrmion phase from 
other phases can be easily determined by the intersects of 
Skyrmion energy line with the other two in Fig. 3. Ap- 
plying dimensional analysis to Eq. (9) or Eq. (11), the 
normalized magnetic field /3 is found to have the same 
dimension as . Therefore one expects both /3ci and /3c2 
to scale with fact also established by numerical cal- 
culation. Restoring proper units, we get the two critical 
fields 
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FIG. 3: (color online) The energies of the three states, i.e., 
(i) helical state, (ii) Skyrmion crystal, and (iii) ferromagnetic 
state. The free energy of ferromagnetic order is set to be 
zero, which is labeled in green. The helical state is energeti- 
cally favored at low magnetic field, while the Skyrmion phase 
emerges at intermediate magnetic field. Ferromagnetism is 
favored at larger fields. 



Bel = 0.2D^/J, Bc2 = 0.8D^/J. (12) 

For Feo.5Coo.5Si thin film, the two critical fields at low 
temperatures were Bd ~ 40mT and Bc2 ~ 80mT^^. Due 
to the co-existence region in the experimental phase di- 
agram, it is not possible to pin down the critical fields 
more precisely. We can still make some estimates based 
on the above analytical formulas of the critical fields. The 
observed spiral wave length of 90nm and the unit cell size 
of a -4.5A gives the ratio D/J = 27r(a/A) w 1/30. The 
critical fields are of order D'^/J = J{D/jf - J/900. 
If we take the observed paramagnetic transition temper- 
ature of ~30K as a measure of J, we would get Be ^ 
30K/900 ~ (1/30)K which in magnetic field unit becomes 
~ (1/30)T, in excellent agreement with the observe field 
ranges of 40-80 mT for the Skyrmion lattice. The criteria 
derived in Eq. (12) is general, applicable to a wide range 
of chiral magnets characterized by both J and D without 
the strong spin anisotropy effects. For such systems the 
Skyrme crystal phase formation is expected in the field 
range of / J . 

The phase diagram sheds light on the properties of 
phase transition as well. Two lines of spiral and SkX 
have different slopes at the critical field /?ci, as shown 
in Fig. 3. Therefore the phase transition between spi- 
ral and SkX phase is first order due to the discontinuity 
of dF/ 8/3 across the phase boundary. For the transition 
from SkX to FM, we carefully examined the derivative 
dF/d(3 of the SkX phase as shown in Fig. 4. At the crit- 
ical field /3c2 = 1.6, dF/d(3 of SkX is clearly nonvanishing. 
However the free energy of the FM phase is set to be zero 
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FIG. 4: (color online) The inter-Skyrmion distance in SkX 
(in red), and the first order derivative of SkX's free energy 
with respect to the magnetic field /3 (in blue). The distance 
increases rapidly near the transition from SkX to FM phase. 
As the derivative is nonvanishing at the critical field /3c2 = 1.6, 
this transition is first order. 



already, so dF/dp = on the FM side. Consequently the 
transition between SkX and ferromagnetic phases is first 
order as well. This is quite consistent with the experi- 
mental results^^. Another interesting issue is the defor- 
mation of Skyrmions at the phase transition from SkX 
to FM. Half of the optimized inter-Skyrmion distance, 
denoted by i?, in SkX is shown in Fig. 4. Deep inside 
the SkX phase, the inter-Skyrmion distance is roughly 
unchanged, which is just the optimal distance 2i?o con- 
trolled by the Dzyaloshinskii-Moriya physics as discussed 
above. However when the magnetic field approaches the 
critical value /3c2, the inter-Skyrmion distance increases 
very rapidly. In this situation, the configuration of the 
individual Skyrmion resembles the one shown in the right 
inset of Fig. 2. The ferromagnetic tail grows significantly 
while leaving the core region almost unchanged around 
radius Rq. The growing FM tail near /3c2 indicates that 
the Skyrmions are separated far away by the intervening 
FM phase, and become more and more dilute when the 
magnetic field approaches /3c2- 

One should be careful to distinguish the SkX with 
R ^ Rq from the FM state. In both states the average 
magnetization is fully saturated, but there is a qualita- 
tive difference in the two phases because of the existence 
of topological defects in the SkX. This explains why the 
SkX— T^FM transition is still first-order (one cannot turn 
off topological defects smoothly), despite the fact that 
the magnetization reaches unity in a continuous manner 
at the SkX/FM phase boundary. A careful examination 
of the Skyrmion size and the inter-Skyrmion distance in 
the thin-film chiral magnet near the upper critical field 
is expected to confirm the field dependence of the two 
length scales Rq and R discussed in this subsection. 
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B. Analogy to Abrikosov Lattice 

In Sec. II, we briefly pointed out the vortexoSkyrmion 
correspondence established by the CP^ mapping (Sec 
Fig. 1). In this subsection, we argue that the corre- 
spondence in fact extends to the case of their respective 
lattice structures. 

Abrikosov put forward the lattice solution of vortices 
in a type-II superconductor of the form^^ 



2J^V — iA — mcr^ z+Bcr2Z = Az. 



(15) 



(13) 



We further assume that the fictitious field produced by 
A is uniform, and choose the Landau gauge = 0, 
Ay = —Hx. The problem of self-consistently deciding the 
vector potential A is reduced to that of a single constant 
H, answering to the first source of non-linearity pointed 
out earlier. As in the typical Landau level problem we 
introduce a plane- wave solution for the y-component, 
z{x,y) = e''^''z(a;), where z{x) obeys 



with some constants Cj . Here ^ is the correlation length, 2J^9a; inax^ z 2J^k+Hx nay^ z+Bcr^z — Az.(16) 



Ix, ly are the inter- vortex separations in the x- and y- 
directions, respectively, with the relation l^ly = h/eB = 
inl"^, ?B=magnctic length. We now pose the question: 
can one construct a spinor solution zskx(a;, y) in the same 
spirit as Abrikosov's vortex lattice solution, such that the 
associated spin configuration nskx = ^skx^'^Skx is the 
Skyrmion lattice? 

We start by writing down the saddle-point equation 
derived by minimizing the total energy, E. The con- 
straint z'''z = 1 is implemented by augmenting the total 
energy with the Lagrange multiplier field A(r): E ^ H = 
E + J (Pr X{z^z — 1). Taking the variational derivative 
5H/ 6z^ = yields the saddle-point equation 

2j(v-iA-iKcryz+2iD{n ■ V)z+{B ■ a)z = \z. (14) 

As a self-consistency check, one can show that the helical 
spin solution (6) indeed satisfies the equation when B = 
0. 

The Skyrme crystal state generates a non-zero 
Skyrmion number that translates into a non-zero effective 
magnetic field, dxAy — dyA^ ^ 0, which in turn should 
generate a Landau-level-type solution. The actual equa- 
tion is non-linear, however, and the underlying Landau 
level structure is not clear. The non-linearity arises from 
two sources. One is that the vector potential A depends 
on the knowledge of the solution z itself. The second is 
the presence of 2iD(n- V)z in the equation. In regard to 
the second issue we are reminded the fact that the SkX 
state exists under a finite magnetic field which partially 
polarize the spins. Numerical calculation shows that the 
mean moment in the SkX state can be a good fraction of 
the full moment^°. It thus appears reasonable to take a 
spatial average of n and obtain n • Vz — > (n) • Vz = 0, 
since the spatial gradient in the two-dimensional lattice 
is orthogonal to the average moment direction along z. 
Such a conclusion will not generally hold for 3D lattice 
or for field tilted away from the .S-direction in the 2D 
lattice. 

Provided the important character of the solution is not 
lost upon the removal of 2iD(n ■ V)z, one can solve in- 
stead of Eq. (14) the following problem 



We have verified that H > corresponds to the 
(over) screening of the external field B by the induced 
field H. 

Introducing the magnetic length Ih = 1/Vh, and 
+ kl'jj, one can derive the solution of Eq. (16) in 
the form zi = (j>n{xk/lH), Z2 = idn(j>n+i{xk/lH) where ^„ 
is the n-th oscillator wave function and d„ is a coeflicient 



dn 



2K^2H{n + l) 



H+B/2J+^{H+B/2J)^ + 8{n+l)K'^H 



■ (17) 



Focusing on the lowest Landau level solution, n = 0, the 
single-particle wave function obtained reads 



^{xk,y) = e 



4>oixk/lH) 

ido (t>iixk/lH) 



(18) 



Wave functions with different fc's are degenerate, and can 
be grouped into a linear combination in the manner of 
Abrikosov solution, Eq. (13): 



/ 2lx \l/2 



E 



pi2TTjy/ly 



-xj/(2i?f) 



.idoV2: 



-</(24) 



• (19) 



Here k is quantized as k = 2-Kj/ly, j=integer, and Xj 

abbreviates x+jlx- Comparison with the vortex lattice 
solution shows that Ih serves as the correlation length 
^, as well as the Skyrmion lattice spacing through the 



condition Ixly 



2-kIIj. The overall constant reflects the 



average normalization (1/Area) / dxdy z^z = (z'I'z) = 1. 

We will now show that the state written down in Eq. 
(19) captures all essential aspects of the Skyrme crystal 
state. Rather than trying to justify the various approxi- 
mations that led to Eq. (19), we regard it as a variational 
state, whose energy can be checked against those of other 
possible spin states. In this way we can construct a phase 
diagram similar to the one shown in the previous subsec- 
tion. 
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For the triangular array of Skyrmions, Cj is chosen 
equal to 1 and i for even and odd integers, respectively, 
and ly equal to \/3 lx/2. Energy per area £ evaluated by 
inserting the variational solution, Eq. (19), into Eq. (5) 
and carrying out the spatial integration reads 



-SkX 



2J 2n^-'-^^ + l±^H 



dl 



dl 



B- 



dl 



-dl- 



(20) 



The extremum condition d£skx/dH = uniquely fixes 
H, hence the magnetic length Ih- One can read off the 
relevant energy scales better by dividing out both sides 
by 4Jk2 = D^/J: 



do 



2V2do 1 l + 3dl 1 ,1 
1 + dg kIh 1 



dl 



dl 



dl' 



\ + 2n^]jb+y/{l + 2nn'ijhY 



(21) 



We used the field in reduced unit b ~ B/{D'^/J) in the 
above expressions. The re-scaling makes it clear that the 
relevant Zeeman energy is indeed D^/J, in agreement 
with the analysis of the previous subsection. 



(b) 



f 

■7/ ; t 
^ \ 1 t 



■ t 1 1 

: \ 1 



, / f r I . . V ^ 



FIG. 5: (color online) (a) A typical Skyrme crystal spin config- 
uration given by Eq. (19) with B — / J , with the optimized 
lattice spacing kIh = •\/3/2. (b) Skyrme crystal spin configu- 
ration obtained by Monte Carlo method from the lattice spin 
model^". 

Figure 5 (a) shows the typical spin configuration given 
out by the CP^ solution zgkx for 6 = 1 (B = 0"^/!). 
Optimizing the energy gives out kIh = ■\/3/2. Spins are 
pointing up, aligned with the B direction, away from the 
Skyrmion center and pointing down at the core. The 
sense of spin swirling (vorticity) in the core region is 
consistent with the right-handed proper screw direction 
of the helical spin phase. Reversing the sign of k, hence 
do — fio in Eq. (19), leads to the left-handed screw and 
a clockwise swirling of spins near the cores. For compari- 
son a typical Skyrme crystal configuration produced from 
the Monte Carlo annealing of the lattice spin model^'^ is 
reproduced in Fig. 5 (b). It is clear that the essential fea- 
tures of the Skyrmion lattice configuration has survived 
the several drastic approximations employed in arriving 
at the Skyrme crystal solution zskx- 



(a) 



(b) 
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FIG. 6: (color online) (a) Self-consistently induced field H 
(in the dimensionless unit k^H) and the ferromagnetic polar- 
ization (n^), against external magnetic field b — B/{D^ /J). 
Taking k,Ih ~ 1/8/2, two different values of u (see text for 
definition), u ^ and u ^ {1/8){D^/J) were used with little 
differences in the results, (b) Energies of helical, Skyrme crys- 
tal, and ferromagnetic spin states against b, both measured 
in units of /J. SkX energy fskx goes up with u while the 
other energies remain insensitive to u. bci and bc2 define the 
two first-order transitions. 



Self-consistently determined H and the average polar- 
ization (n^) = (zV^z) = (1 — c?o)/(l + ^0) against the 
external field b are shown in Fig. 6 (a). It is seen that H 
tends to zero as the external field drives full polarization 
of spins, (riz) ~ 1. The core size and the Skyrmion spac- 
ing, both of which are fixed by Ih, diverges accordingly. 
In reality the divergence is cut off by a first-order transi- 
tion to the energetically more favorable FM state at the 
upper critical field Bc2- 

Shown in Fig. 6 (b) are the energies of the three 
main competing spin configurations - helical, Skyrme 
crystal, and ferromagnetic spin states - plotted against 
b. The SkX state gives the lowest energy regardless of 
the field strength, in contrast to both Monte Carlo^° and 
experimental^^ findings showing the ground state evolu- 
tion helical— >SkX^FM with increasing B. We also find 
that the energy is insensitive to the crystal structure of 
the Skyrmions being a square or a triangular lattice. 

Both these problems can be remedied by re-visiting 
the constraint z^z = 1, which is so far imposed at the 
crudest, mean-field level. Helical and FM spin solutions 
obey the constraint exactly anyway, but the SkX solution 
does not. An inclusion of the potential term u(z^z — 1)^ 



8 



{u > 0) in the free energy density (5), which imposes the 
constraint at the local level, would increase the energy 
density of the SkX solution by u[((zg^j^zskx)^) — !]• As 
seen in Fig. 6 (b), a suitable choice of u restores the cor- 
rect sequence of ground states punctuated by two critical 
fields, Bel for hehcal— >-SkX and Bc2 for SkX— >-FM transi- 
tions. According to Fig. 6 (a), H remains close to unity 
in the whole SkX region, implying the inter-Skyrmionic 
spacing Ih comparable to the period of the spin spiral 
A ~ K~^. This is indeed the case both in the Monte 
Carlo calculation^^ and experimentally^'^^. The chance 
of a larger Skyrmion spacing as H becomes very small is 
pre-emptied by the transition to a ferromagnetic phase 
with a lower energy at B = Bc2- Finally, the average 
((zgj^-j^zskx)^) calculated at a number of {B,H,k) values 
is in favor of the triangular lattice having a lower energy 
over the square lattice. This is the same criterion used 
by Abrikosov^^ and subsequent workers^" in identifying 
the most stable lattice structure of vortices by calculat- 
ing the average of \tp\'^ in the U(l) GL theory. Unlike in 
Abrikosov's case, the average ((zgj^j^zskx)^) depends on 
the {B, H, k) values used and are not universal. 
With zskx at hand, the induced field distribution 

-{2TT/ly){x/Qd^{zl^^^zsi,x) (22) 

can be worked out^^. Figure 7 (a) gives the field distri- 
bution, along with the distribution of the z-component of 
the local magnetization n^{x,y) in Fig. 7 (b). The field 
intensity reaches a maximum at the Skyrmion core and 
nearly equals zero in the FM background. The spatial 
average of Hs\^x{x, y) is approximately equal to H used 
as an input in zskx- 




FIG. 7: (color online) (a) Induced field distribution 
Hskxix,y) from the Skyrmion lattice solution zskx- Bright 
area is the Skyrmion core where the intensity is the maxi- 
mum, (b) Distribution of n^{x,y) = Zgj^-j^cr^zskx- Dark area 
means reversed spins at the Skyrmion core. 

We should emphasize that in the 3D chiral magnets it is 
the conical phase, with the propagation vector along the 
z-direction, having the lower energy and hence occupying 
the lower temperature side of the phase diagram instead 
of the SkX state, under moderate magnetic field^"*. The 
conical phase cannot exist in a 2D system, and this prob- 
ably contributes to the SkX state being realized at low 



temperature for thin-film samples. It still remains an in- 
teresting challenge how one can understand the crossover 
in the behavior from 2D to 3D systems. 

C. Comparison to earlier work 

Bogdanov and collaborators have extensively investi- 
gated the possibility of the Skyrmion formation in chi- 
ral magnets on the basis of Ginzbug-Landau models 
with Dzyaloshinskii-Moriya interaction^^'^'^. In partic- 
ular the field dependence of the phase diagram presented 
in Ref. 13 (Fig. 9 in their paper) correctly captures the 
spiral— >SkX-^-FM phase change recently investigated by 
Monte Garlo method^°, experimentally^^, and in the pre- 
vious two subsections. The analysis of subsection III A 
amounts to their "circular-cell" approximation. 

In Ref. 13 the nature of the phase transition from 
SkX to FM was not determined conclusively. Based on 
our analysis presented in Sees. Ill A and IIIB, it is 
clearly first-order. Physically this should be clear since 
there is no way to turn off topological defects in a con- 
tinuous fashion. Nor is it possible to fuse a Skyrmion 
with an anti-Skyrmion to annihilate them, since the sys- 
tem under consideration consists only of one species of 
Skyrmions^^^'^*^'^^. Experiments also support the first- 
order phase boundary 

In their analysis, Bogdanov et al. assumed a two- 
dimensional structure homogeneously extended along the 
third direction. Our model is explicitly two-dimensional. 
With the recent experimental input^ we know that for 
known 3D chiral magnets such as MnSi, in the field range 
where the SkX phase (a.k.a. A-phase) is found, the low- 
temperature state has the conical spin structure. In order 
to stabilize the Skyrmionic phase it is therefore essential 
to suppress the three-dimensionality as well as to turn 
on a magnetic field. The present investigation fills in the 
gap that existed between the early theoretical work and 
the new insights offered by recent experimental progress. 

Finally, our paper presents a new framework for deal- 
ing with chiral magnetic systems. The CP^ formulation 
allows us to understand the Skyrmion lattice formation 
in direct analogy to the Abrikosov vortex problem in su- 
perconductors, and to the Landau-level physics. It is 
our belief that the new formulation can be used in con- 
junction with the conventional classical vector theory in 
studying other aspects of chiral magnetism such as their 
dynamics. 

IV. DISCUSSION 

Although the detailed examination of the thermal ef- 
fects at nonzero temperature is outside the scope of this 
paper, we can make qualitative assessment about the 
thermal fiuctuations for the Skyrme crystal state. Both 
translational and spin rotational symmetries are bro- 
ken in the SkX state, and one accordingly expects two 
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types of Goldstone modes: phonon-like displacements of 
Skyrmion positions and spin waves. 

Applying the Lindcmann criterion, the melting of the 
Skyrmion lattice will take place if the thermal fluctua- 
tion in the center-of-mass position of a Skyrmion becomes 
comparable to the inter-Skyrmion spacing. As discussed 
in the previous section, the typical magnetic energy dif- 
ference per spin for the three magnetic phases - helical, 
SkX, and FM - arc of the order D"^ / J . A shift of the 
Skyrmion position by one lattice unit increases the en- 
ergy per spin by a typical amount D'^/J, or JZ^^ if the 
linear size of the Skyrmion Ih ~ ^ J/D is used. 
The energy increase associated with the Skyrmion dis- 
placement by X should then be J x (x/Ih)^- Applying 
the equipartition theorem we obtain J x {x'^)/Ih ~ T at 
a temperature T, an d the melting temperature ~ J 
from the condition (x^) ~ In- 

Another possible disordering mechanism is the spin 
fluctuation. In the ground state each Skyrmion unit cell 
exhibits identical spin configuration. If the spin orienta- 
tion of one particular Skyrmion unit cell is completely re- 
versed, the cost in Zeeman energy is roughly B times the 
number of spins, or Blfj. Noting that B in the relevant 
SkX phase is /J, we conclude that the single- Skyrmion 
spin flip has the energetic cost of J. From this we con- 
clude that the coupling energy of adjacent Skyrmions 
must be J X (A^/tt)^ for a small angle difference 
of the two nearby Skyrmion spin orientations. By in- 
voking equipartition theorem again, we arrive at a spin- 
melting temperature ~ J. In conclusion, destroying 
the global ordering of the SkX state by either center-of- 
mass displacement or spin fluctuation requires a temper- 
ature of order J, even though the magnetic field energy 
difference between the various phases is only D'^/J. It 
is the large size of the Skyrmion, of order (J/D)^, which 
compensates for the small magnetic energy scale and ren- 
ders a sizable melting temperature. 

To summarize, we have provided a theoretical treat- 



ment of the Skyrme crystal phase recently observed in 
the thin film of chiral magnet Fco.sCoo.sSi^^. Two in- 
dependent constructions of the Skyrmion lattice state 
were used to calculate the energy of this state and 
compare it against those of other competing states, 
namely helical spin and ferromagnetic spin, as the per- 
pendicular magnetic field strength is increased. The 
resulting phase diagram giving the successive ground 
state evolution of helical^SkX^-FM phases is in ac- 
cord with the Monte Carlo simulation^° and the thin-film 
experiments^. We gave a general argument for the two 
critical field strengths separating the SkX state from the 
helical state on the low-field side and from the FM state 
on the high-field side to be of order D^/J. Other spi- 
ral magnets with similar spin structures are expected to 
exhibit Skyrme crystal phase for field ranges of D'^/J if 
they are made in the thin-film form. We also gave a qual- 
itative argument why the melting of the Skyrme crystal 
should occur at T„ ~ J despite the small magnetic en- 
ergy scale D"^ /J needed to stabilize the SkX phase. As an 
interesting by-product, we derived the variational wave 
function corresponding to the Skyrme crystal state in the 
form. 
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